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I. INTRODUCTION 

Quantum entanglement is a key resource in the 
theory of quantum information, including quantum 
cryptograph^, quantum communication^ and quantum 
computations.—. In the last few years significant efforts 
have been undertaken in order to engineer devices capa- 
ble of generating maximally entangled states. In partic- 
ular, in view of their integrability and scalability, consid- 
erable attention has been devoted to the generation of 
entanglement for electrons in solid state systems (for a 
review see Refs. HH). Proposals aim at producing pairs 
of entangled electrons either in the spin or in the or- 
bital degrees of freedom. Some of them are based on the 
generation of Bell states by means of electron-electron 
interaction^^ ! 10 ! 11 ! 12 . Only recently it has been shown 
that in a mesoscopic multiterminal conductor entangle- 
ment can be also produced in the absence of electron- 
electron interactio n 13 : 14 ' 15 ! 16 ) 17 : 18 : 19 and through time- 
dependent potentials' 2 *^ 2 -. 

Unfortunately, electronic entangled states are fragile. 
Indeed, it is well known that the unwanted interaction 
with the environment might lead to loss of quantum co- 
herence and eventually to the destruction of the entan- 
glement between the electronic quantum states. In view 
of experimental implementations, it is therefore of cru- 
cial importance the search for simple procedures able to 
assess the entanglement actually generated in a given de- 
vice (the entangler) under non vanishing environmental 
couplings. 

Several strategies have been proposed in order to de- 
tect electronic entanglement in mesoscopic conductors. 
One of these schemes makes use of an additional beam 
splitter, to be attached to the entangler generating pairs 
of spatially separated entangled electrons, which reveals 
the presence of entanglement through a measurement of 
the correlations of the current fluctuations in the exit 
leadsS 3 -. The possibility of evaluating a lower bound for 
the entanglement of formation was later demonstrated in 
Ref. Inland, recently, a generalization in terms of a Hong- 
Ou-Mandel interferometer was presented 2 ^ The influ- 
ence of Rashba spin-orbit couplin g 26 ! 27 , rotating mag- 



netic fields 2 ^ and dephasing contacts 2 ^ were also inves- 
tigated. In addition, it has been shown that signa- 
tures of entanglement in this prototype setup can be 
detected through the full counting statistics 3 ^. A sec- 
ond scheme makes use of Bell inequality tests, which 
have been formulated both in terms of current corre- 
lations at different tcrminalsi^ii 3 -^^^^^ and by 
resorting to the full counting statistics 1 ^ 2 -. More re- 
cently, the following other schemes have been consid- 
ered. For a superconductor-dot entangler, spin current 
correlations have been shown to give a characterization 
of the entangler operation^ 3 -. A one-to-one relation be- 
tween the entanglement production rate of electron-hole 
pairs and the spin-resolved shot noise has been derived 
in Ref. [H] in the case where the elastic scattering is 
spin-independent. A proposal for detecting energy en- 
tanglement in a normal-superconducting junction, mak- 
ing use of a Mach-Zehnder interferometer, has been also 
presented 3 -^. Furthermore, a quantum state tomography 
scheme has been put forward 3 ^ 

In this paper we shall show that it is possible to im- 
plement the idea of entanglement witness operators 37 to 
detect entanglement for electrons in multiterminal meso- 
scopic conductors afflicted by "noise". The concept of 
entanglement witness has been introduced only recently 
in quantum information theory^ and experimentally ap- 
plied to optical systems 3 ^. It has been shown that, given 
a class of states (containing both separable and entangled 
states), it is possible to construct an Hermitian operator, 
called witness, having the remarkable property that its 
expectation value is positive for all the separable states 
contained in that class. If this is true only for separable 
states, the witness is said to be optimal. As a conse- 
quence, all entangled states are singled out by negative 
expectation values of the optimal witness operator. Al- 
though the entanglement witness detection is not a mea- 
sure of entanglement, it represents an efficient method to 
detect entangled states. 

In order to implement optimal witness detection for 
electrons in mesoscopic structures one needs: (i) to iden- 
tify the class of states that is generated by the given 
setup, (ii) to construct the optimal witness operator, (iii) 
to be able to measure its expectation value. Given a 
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physical setup and a proper phcnomcnological model for 
the "noise" afflicting the device, the class of states can 
be easily identified. Also the construction of the optimal 
entanglement witness does not present much difficulties 
once the separability criterio n 41 ' 42 can be used. Unfor- 
tunately, the measurement of the witness operator is not 
always an easy task. Furthermore, the witness operator 
can be measured only when it does not depend on the un- 
known parameters describing the presence of "noise" . In 
that case we shall demonstrate, for different setups, that 
the mean value of the witness operator can be determined 
through the measurement of current cross-correlations. 

The effectiveness of the entanglement witness stems 
from the fact that the set of cross-correlated measure- 
ments required to determine its expectation value is fixed, 
irrespective of the degree of "noise" present, and the 
number of such measurements is small. Typically, from 
3 to 5 analyzer's settings are necessary, corresponding to 
a number of current cross-correlations of the order of 10. 
In this respect it might be preferable to a Bell inequality 
test, which, in the presence of "noise" , demands an un- 
known and possibly large number of analyzer's settings 
to be tested in order to search for violations. Further- 
more, the Bell inequality test is only a sufficient criterion 
for the detection of entanglement 4 ^, while the witness, if 
optimal, is both sufficient and necessary. Witness entan- 
glement detection might also be preferable to quantum 
state tomography, especially when one is not interested 
in acquiring information over the whole density matrix of 
the system, but simply wishes to know if some entangle- 
ment has survived at the exit of the setup. 

The paper is organized as follows. In Sec. [TT] we in- 
troduce the concept of witness operator and discuss how 
the latter can be easily constructed by using the separa- 
bility criterio n 41 ' 42 . In Sec. IIIII wc describe the details of 
witness entanglement detection by applying it to a proto- 
type setup where electron spin singlet and triplet states 
are initially generated and subsequently corrupted by de- 
coherence due to spin dephasing and spin relaxation in 
the wires. We then discuss under which conditions the 
measurement of the witness operator can be achieved 
and we show that the expectation value of the witness 
operator can be expressed in terms of current cross- 
correlators. We also explain why entanglement witness 
detection might be more efficient than usual Bell inequal- 
ity tests. In Sec. [TV] we then proceed to construct the wit- 
ness operators for two realistic setups generating particle- 
hole entanglement. Specifically, in Sec. II V Al we consider 
the quantum Hall bar system proposed in Ref . |I3| while in 
Sec. IIVBI we discuss the electronic Hanbury-Brown and 
Twiss interferometer proposed in Ref. HH Sec.|V]is finally 
devoted to the conclusions and future perspectives. 



II. THE WITNESS 

Let p be the state of a bipartite system acting on 
the Hilbert space Tt = TCa ® Hb and let us denote with 



A and B the two subsystems. The state p is called 
entangled^ (or non separable) if it cannot be written 
as a convex combination of product states, i.e. as: 

p = ^2pk\ak,b k )(a k ,b k \ (1) 

k 

where pk > and \ak,bk) = \ak)A <8> \bk)B are product 
vectors. Conversely, the state p is called separable (or 
not entangled) if it can be written in the form given in 
Eq. 0. 

It is well known that, limited to the case of TC = C 2 ®C 2 
or H = C 2 <£> C 3 , it is possible to determine whether p is 
entangled or not by calculating the partial transpose (i.e. 
the transpose with respect to one of the subsystems). 
In fact, according to the separability criterio n 41 ' 42 , the 
state p is separable if and only if its partial transpose 
is positive. Alternatively, one can distinguish between 
an entangled state p and all the separable states by in- 
troducing the concept of entanglement witness^. Given 
an entangled state p, an entanglement witness is indeed 
an Hcrmitian operators W having the remarkably propri- 
ety that Tr(Wp) < and Tr(Wp scp ) > for all separable 
states Psep- The negative expectation value of W is hence 
a signature of entanglement and the entangled state p is 
said to be detected by the witness. 

Given a class of states p\ (containing both separable 
and entangled states) we define as optimal^ the witness 
operator W op t that is able to detect all the entangled 
states present in the class pa- As a result, whenever 
a class of states p\ and the optimal entanglement wit- 
ness for that class are given, it is possible to determine 
whether a state p £ p\ is entangled by simply checking 
if (W^pt) = Tr(Woptp) is negative. 

As one can expect, finding an optimal witness for a 
given class of states is, in general, not a trivial task. 
However, the optimal witness construction can be greatly 
simplified by resorting to the separability criterion 4 ^. In 
that case, one finds that the optimal witness operator 
reads: 

Wopt = |0)(0| TB (2) 

where \4>) is the normalized eigenvector corresponding to 
the negative eigenvalue of the partial transpose of the 
density matrix, i.e. p^J 3 (Tb represents the partial trans- 
position with respect to the subsystem B). 

In view of an experimental measurement the mean ex- 
pectation value of the optimal witness (Wopt), it is nec- 
essary to decompose the optimal witness into operators 
which can be measured locally. In this paper, for exam- 
ple, we shall use the following decomposition: 

Wopt = ^2 K nmOVi®C m , (3) 
n,m— 0,x,y ,z 

where o~q is the unity 2 x 2-matrix 1 and a n , with 
n = x,y,z, are the Pauli matrices. Notice that, accord- 
ing to this procedure, the average value of the optimal 
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FIG. 1: Prototype setup consisting of an entangler (yellow 
region) connected to terminals A and B through two con- 
ductors (green regions) characterized by transmission ampli- 
tudes tp^€'' i ' Aa and tBe l ^ Ba . The current is measured spin- 
selectively in A and B along spin quantization axes repre- 
sented by the vectors a and b. 



witness operator can be calculated only if the coefficients 
Knm are known. Moreover, it is highly desirable that 
the decomposition given in Eq. (|3|) contains a minimal 
number of terms. Since we shall only consider optimal 
witness operators, in the following the word optimal will 
be understood. 



III. PROTOTYPE SETUP 

For defmitcness let us consider first the prototype sys- 
tem presented in Fig. ([I]) composed of an entangler con- 
nected, through two conductors (green wires in the fig- 
ure), to two outgoing terminals, A and B, where electric 
current can be measured spin-selectively along specified 
quantization axis (a and b for upper and lower electrodes, 
respectively). The entangler is a device which produces 
maximally entangled (Bell) states: 



I*) = 



n 



|4 T (£)4u(£) + (-iya f M (E)a^(E) 



0<E<eV 



V2 



|0>, 
(4) 

where s = 1 (0) for spin singlet (triplet) and a? Aa (E) 

( a Bo(E)) is a creation operator for an electron with spin 
a in the upper (lower) conductor at energy E. The prod- 
uct is taken over energies ranging from the chemical po- 
tential of the outgoing terminals (set to zero) and the 
energy eV, while the state |0) represents the filled Fermi 
sea. We denote by t\ (fa) the transmission amplitude 
of the upper (lower) conductor, assuming it as energy- 
independent in the range < E < eV. Furthermore, we 
shall include the presence of "noise" in the conductors 
through two different models. In the first one (Sec. lIII A")) 
we only assume dephasing of the spin degree of freedom, 
while in the second (Sec. IIII Bp we also introduce spin- 
relaxation. 



A. Random phases accumulated in the wires 

We first consider a rather simplistic model where the 
"noise" in the transmission wires is simply an acquired 
random phase. Specifically, we introduce spin-dependent 
random phases (</>Aa and 4>Ba), so that the overall ef- 
fect of the conductors on the electronic operators can be 
described by the transformation: 



al(E)-,Ue^al(E), 



(5) 



where i — A, B and a =|, J,. The density matrix of the 
outgoing state, i.e. of the electrons arriving at the exit 
terminals, is given by 



1 



u(0)eV 



(0) 



0<E<eV 



with 



P{E) = - 



where 2 a 
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(7) 

6ax, £b = 0B| — </>bj. & n d v{E) is the 
density of states, which we assume to be independent of 
energy in the range considered. The subscript E indi- 
cates that the matrix refers to the subspace relative to 
energy E. Note that, as it should be, entanglement is 
not affected if up and down electrons in the same con- 
ductor accumulate an equal random phase, since in that 
case xa = xb = 0. We assume a Gaussian distribution 

for the random phases, i.e. P{4>i a ) = j^-nu e w ^ n 
uj being the width of the distribution, termed dephasing 
parameter below. The phase-averaged density matrix (p) 
is given by replacing, in Eq. ([7|), e ±l ( XA ~ XBS> with e -2 " . 

In this setup the two subsystems are identified with 
conductors A and B, so that the partial transpose with 
respect to B is given by (p TB )ij-ki = (p)u-kj, where i and 
k arc spin indexes for A, while I and j are spin indexes 
for B. In order to construct the witness we perform the 
partial transpose with respect to the subsystem B on 
the phase-averaged density matrix and find that the only 



negative eigenvalue is A 



, with eigenvector 



0<E<eV 



n ^[-4 T (^)4 T (^)+ 



4n 



(E)}\0). 
(8) 

The fact that A is always negative proves that the outgo- 
ing state is entangled for every value of the dephasing pa- 
rameter lu smaller than infinity, irrespective of the value 
of the transmission amplitudes tx and t^. The witness 
is then given by W = ® 0<E<eV We, where 



W E = - 



10 

o (-i) s+1 o 

(-f) s+1 

Of 



(9) 
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so that, by construction, the mean value of the witness 
operator (W) = Tr(Wp) equals A 

It is important to stress that, in order for the witness to 
be of use, W has to be independent of unknown param- 
eters like. The decomposition, containing the minimal 
number of terms, for the witness operator in Eq. ([9]) is 
given by^ 



W E = + 



1 E 

k—x,y 



Ok 65 o k 



(10) 



As it is evident from Eq. (fTUJ) , operatively the determi- 
nation of the mean value of W requires the measurement 
of the spin correlators (a - o ®b ■ o) for a few sets of a 



J 



and b. In the following, we shall see that under certain 
conditions the average value of such correlators can be 
given in terms of current cross-correlators. 



To this aim, let us consider the correlator 

(Q Aaa{tm)Q Bbpitm))^; where Q Aaa{tm) \Q sg^^m) 

is the number of particles detected over a measuring 
time t m in A (B) with spin a (/?) along the direction a 
(b). Intuitively, such a correlator is proportional to a 
coincidence counting, at least, in the situation where no 
more than one particle per port arrives at the measuring 
apparatus within the measuring time t^^. Therefore, 
one expects the following relation to hold^ i n ' 13 i 31 ' 49 : 



(a ■ o 



[QAat(t m ) - QAal(im)][Q B £ T (im) — Q B ^(*m)]} 
[Q A al(t m ) + Q A Sl(tm)][Q B ^(t m ) + Q B - bi {t m )})' 

I 



(11) 



with correlators defined as^ 

(Q A a*(t m )Q BSf3 (t m )) = P dtj m dt'{I S Aa {t)I%{t')). 

(12) 

In Eq. (TT2"]) I% a (t) is the current operator, defined in Ap- 
pendix [XI] for electrons arriving at terminal A with spin 
a along the direction a, and analogously for I b B3 {t'). We 
now make use of the fact that the correlator in Eq. JT5|) 
can be expressed in terms of the zero- frequency current- 
correlator defined as 

/oo 
dt(SI s Aa (t)5I b B0 (O)), (13) 
-OO 

with Slf a (t) = lf a (t) - (lf a ). The actual expression de- 
pends on whether the measuring time t m is larger or 
smaller than the correlation time r c = h/eV, which 
represents the time spread of the electron wave packet. 
Namely, for t m ^> r c one gets^ 

(QA Sa (t m )Q B ^(t m )) ~ t 2 m (li a )(I Bf3 ) + t m s%(0) (14) 

and, for f m C t Ci one gets 

(Q A aa(tm)Q BSf3 (t m )) * t 2 m [(I Aa )(ll p ) + ^4(0)]. 

(15) 

At zero temperature and for short measuring times 
(t m <C t c ), the spin correlator for our setup turns out 
to be£i 



(a-a<E)b-(f) 



eV3$(0) + sg(0) - sg(0) - sf»(0) 
2h (I%+I S Ai )(I%+lli) 



(16) 



By using Eq. 
witness is 



TOl. we find that the mean value of the 



(W) 



"2 6 



(17) 



which equals the values found for A. On the contrary, in 
the limit of large measuring times (t m 3> t c ) we obtain: 



la ■ a ®b ■ o) = 



1 ^Tt(°) + - s fi(°) ~ s 4t(°) 

2<m (Jlt + J l)(4t+4i> 



(18) 

which is always vanishing, even in the tunneling regime. 

This is due to the fact that, in our setup, current cross- 
correlations of Eq. (fl~3"|) are inherently a result of two par- 
ticle scattering processes. More precisely, current corre- 
lations are proportional to | i a | 2 1 *b | 2 likewise the product 
(I Aa )(I B p)- In contrast, in fiefs. [Tllfl3l . for example, the 
cross-correlations result from partition noise and, in the 
tunneling regime, are much larger than the product of the 
currents (times < m ) provided that the measuring time is 
smaller than the average time spacing between tunneling 
events. It results that for large measuring times one finds 
(W) ~ j, which would be the result obtained averaging 
over a fully mixed state. We conclude that the witness 
operator can be measured only for measuring times t m 
much smaller than the correlation time r c . Nevertheless, 
entanglement can be detected even beyond the latter con- 
dition. Indeed, for arbitrary measuring times t m the spin 
correlator can be written as: 



(a-o®b-o) = a(t m ,T c ) 



4*(0) + S |J(0)- S fJ(0)- S f*(0) 



(19) 
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where a is, in general, a function of both r c and t n 
following simple relation can then be derived: 



(W) = 2a(t r 



= )TcA + - [1 



2a(t m ,T c )r c 



The 



(20) 



It turns out that, for 2ar c > 1, (W) < implies A < 
so that entanglement can be witnessed as long as (W) is 
negative. Furthermore, since the spin correlator in Eq. 
(fT9|) depends explicitly on t m , its measurement requires 
time-resolved detection, i. e. energy filtering. This is not 
the case in the limit t m <C t c . 

It is interesting to compare the capability of detec- 
tion of entanglement provided by the witness with the 
one relative to a Bell inequality test. In the following 
we choose the well known Clauser-Horne-Shimony-Holt 
(CHSH) inequality^±ML^: 



£ 



E(a,b) + E(a' ,b) + E(a' ,b) - E(a' ,b') < 2, (21) 



where E(a, b) = (a-<j®b-<j) is the spin correlator given in 
Eq. p6[) and a,a',b,b' denote arbitrary directions along 
which the spin is measured. In the absence of dephas- 
ing, for the case of the Bell state given in Eq.(|D), the 
CHSH inequality is maximally violated by choosing the 
following set of directions: 

a = (1,0,0) 6=-L(i,0,l) 

a = (0,0,1) b' = -^(1,0,-1). (22) 

Indeed one finds that £ = £ max = 2 y/2 and, as a result, 
the entanglement of the Bell state is detected. Notice 
that, in order to experimentally test the inequality, for 
every set of directions one needs to measure 4 spin corre- 
lators and, for each of them, 4 current cross-correlators. 

Let us now assume the presence of dephasing in the 
wires. It is easy to verify that, for t m <C r c , £ = V2(l + 



e 2u ). As a result, for value of the dephasing parameter 
uj > J i| ln(\/2 — 1)| the Bell inequality is not violated 
even thought, according to the previous witness detection 
analysis, the state exiting the leads is still entangled. 

It has to be noticed, however, that if the dephasing 
parameter uj is known, it is possible to find the proper set 
of directions leading to violation of the CHSH inequality. 
In practice, since the value of uj is unknown, in order to 
find a violation of the inequality, one needs to try several 
measurements of current cross-correlators over different 
sets of directions. On the contrary, by measuring the 
witness given in Eq. (|10p one needs to perform only 12 
measurements of current cross-correlators, 4 for each spin 
correlator, with fixed configurations irrespective of the 
degree of dephasing in the conductors. 

As a final remark, we notice that the measurement of 
the mean value of the witness gives information on the 
degree of the dephasing affecting the wires. 

B. Spin relaxation and dephasing 

We shall here analyze the same system considered in 
the previous section, now assuming that during the prop- 
agation in the wires the electron spins are subject to re- 
laxation in addition to dephasing. The former might be 
caused by the presence of magnetic impurities, nuclear 
spins, or the spin-orbit coupling in the wires^. We shall 
adopt a phenomenological model of decoherence along 
the lines of Ref. [23|, by resorting to the Bloch theory^ 
and by expressing the density matrix at each exit leads 
as: 



Pi(t L ) = Ai(t L )\pi(0)] 



(23) 



with i = A, B. Here = L/v-p is the ballistic transmis- 
sion time (L is the length of the wires and vp is the Fermi 
velocity) and we have introduced the superoperator: 



J 



Ai(tz)[/>] = 



Ptt (l + e~W Tl ) + PU (l-e"W Tl ) 



2e -WT 2 



Pit 



2e ~tL/T 2 



Pn 



ptt (i- e -*- /Tl )+m(i + e- ti/Tl ) 



(24) 



where T\ is the relaxation time and T2 is the dephasing 
time. By assuming that the decoherence and relaxation 
processes act on each wire independently, we can readily 
write the density matrix of the outgoing state as follows: 



p{t L ) = {A A {t L )®A B {t L ))[p{0)} 



(25) 



where p(0) is the density matrix of the initial entangled 
state given in Eq. It is a matter of simple calcula- 



tions to show that 
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where d± — (l ± e _2tl -/ Tl ). The witness operator can 
now be determined by applying the method described in 
Scc. lIII Al We find that it is equal to the one reported in 



FIG. 2: (Color online) Mean value of the witness as a function 
of inverse relaxation (th/Ti) and dephasing (£l/T2) times. 
The gray plane corresponds to (W) — 0. A typical GaAs 
structure has L ~ lfim and vf ~ 10 4 — W'm/s, so that we 
obtain tn w 10 — lOOps. Typical decoherence time found in 
GaAs are of the order T2 ~ 100ns — 1/is. It follows that 
t L /T 2 « 10~ 5 - 10~ :i . Typically Tj > T 2 . 

Eq. (JHJ) and its mean value reads: 

(W) = Tr (Wp(t L )) = - (l - e~ 2 ^ - 2e~ 2 %) . (26) 

In Fig. [2] wc report a three-dimensional plot of (W) 
as a function of tz/Ti and £l/72. To highlight the 
region of the surface which corresponds to entangle- 
ment, we have plotted a gray plane in correspondence to 
(W) = 0. Entanglement is present when the orange sur- 
face is below the gray plane. For large relaxation times 
(ti/Ti — > 0) we find that (W) is always negative: in 
agreement with the results of Sec. IIII A[ we find that 
for arbitrary degrees of dephasing there is always en- 
tanglement in the outgoing states. However, for finite 
relaxation time, there always exists a threshold for the 
dephasing above which entanglement is lost during the 
transmission in the wires. Indeed it is easy to verify that 
for t L /T 2 > || In i (l - e - 2tL / Tl ) I the mean value of the 
witness becomes positive. Note furthermore that relax- 
ation processes have a much less disrupting effect with 
respect to dephasing processes, due to the prefactor 2 in 
Eq. (EH). 

Let us again compare the performance of the witness 
to the Bell inequality test. We find that, choosing the 
set of directions given in Eq. (|22p . the CHSH parameter 
reads 

£ = |V2 (V 2 ^T + e~ 2 %) |. (27) 

In Fig. [3] we plot the CHSH parameter as a function of 
t^/Ti and thjT%. We draw a gray plane in correspon- 
dence to £ = 2 in order to highlight the region in the pa- 
rameter space (above the plane) for which entanglement 



FIG. 3: (Color online) CHSH parameter £ as a function of 
inverse relaxation (fi/Ti) and dephasing (t^/Tb) times. The 
gray plane corresponds to £ — 2. 

is detected. By comparison with Fig. [21 it is evident that 
there is a large region in the parameter space for which 
Bell inequality fails to detect entanglement. Notice that 
also in this case £ can be maximized by exploring differ- 
ent sets of angles. 

IV. PARTICLE-HOLE ENTANGLEMENT 

A. Edge states in a quantum Hall bar with a 
quantum point contact 

It is interesting to apply the method of the witness to 
some realistic systems. Let us now consider the particle- 
hole entanglement of edge channels in the quantum Hall 
bar system proposed in Ref. [H. The structure, depicted 
in the inset of Fig. 01 consists of a quantum Hall bar di- 
vided in two halves by a quantum point contact (QPC). 
A large magnetic field is applied perpendicular to the 
plane, so that two edge states are allowed. Electrons 
are injected from a contact, represented by a red rectan- 
gle, biased at a voltage V. They propagate to the right 
through the two edge channels (labeled by 1 and 2), and 
are transmitted to the right-hand-side with amplitudes 
matrix t and reflected in the left-hand-side with ampli- 
tude matrix r. In the tunneling limit, the outgoing state 
can be written as^ 

I*) = II (v^l*> + VT^filO}) , (28) 

E 

where 

m = n- 1/2 k{ibl\d) (29) 

is an entangled state expressed in terms of particle oper- 
ators 6p = {b pl ,bp 2 ) for the right-hand-side and hole op- 
erators frjj = ^2) f° r th e left-hand-side. In Eqs. I|28p 



7 



and (|29l) . = Tr^y-y', 7 = a y ra y t T and |0) represents a 
new vacuum state consisting of particle-filled edge states 
in the left-hand-side. Two additional gate electrodes, 
which provide a local mixing of the channels, are placed 
on the left and on the right with scattering matrix [7l and 
Un, respectively. Current cross-correlators, which deter- 
mine the spin correlator through Eq. pip , are measured 
between the exit ports L and R. 

Let us consider the following parametrization of 7^: 



7 = e ie y/T(l-T) 



gi(<^L2+0Bl 
_gi(0Ll+0R 



cos£ 
> sin£ 



o i(0L2 + 0R2) 
a i(0Ll+0R2) 



sin£ 
cos£ 
(30) 

which implies identical transmission eigenvalues T for 
the two channels. The possible presence of dephas- 
ing processes is accounted for by introducing Gaussian- 
distributed random phases 0l(r)j with equal width to. £ 
is the channel-mixing parameter of the QPC. We find one 
negative eigenvalue 



_ V2(l - 7] 2 ) - y/1 + 6r? 2 + if + (-q 2 - l) 2 cos 4^ 

(31) 

for the partially transposed phase-averaged density ma- 

2 

trix, where 77 = , showing that entanglement exists 
only within some ranges of the parameters £ and r/. It 
is interesting to note that A is related to the concurrence 
C calculated in Ref. H3, namely A = — C/2 for positive 
values of C. 

The witness operator, constructed with the separabil- 
ity criterion as in Sec. IIII1 in this case depends on such 
parameters and cannot be used to detect entanglement 
in the general case. If, however, the mixing parameter £ 
is known or controllable, then a witness independent of 
parameters can be constructed for maximum (£ = 7r/4) 
or no mixing (£ = 0,7r/2). In this latter case We takes 
the form of Eq. (|10[) , while in the former we find 



W E 



1 



{1 <8> 1 



a x + a x <g> <r z } . (32) 



Even for this system the possibility of measuring the 
witness depends on the relative position of the time scales 
involved: measuring time i m , correlation time r c = h/eV 
and tunneling time r tun = h/(2eVT). Namely, both 
in the maximum and minimum mixing conditions, for 
t m ^> t c we find that the mean value (W) is measurable 
and equals the value of A given in Eq. (|3ip only in the 
tunneling regime (t m <C rtun). Nonetheless, for t m <C r c 
the witness can be measured irrespective of the value of 
T . Specifically, in the case of no mixing, we find that 
{W) = —r] 2 /2, showing that, irrespective of the degrees 
of the dephasing, the state is always entangled. Interest- 
ingly, this result equals the witness mean value obtained 
in Sec. IIII Al On the contrary, in the case of maximum 
mixing we find that (W) = 1/4(1 — 2i] — if). We plot 
this result in Fig. |4] as a function of 77. Notice that the 
state is entangled only for values of rj > v3 — 2\/2 ~ 0.4. 
This result can be compared with the prediction given by 



a Bell inequality test that assures that entanglement is 
present only for r\ > 1/V2 ~ 0.7 (see the discussion in 
Ref. I&il) . Unfortunately, no parameter- free witness could 
be found for arbitrary mixing £. 




0.25 



-0.25 



-0.5 



FIG. 4: Witness mean value as a function of the dephasing 
parameter r\ for particle-hole entanglement. Inset: the system 
consists of a Hall bar divided in two halves by a QPC. A large 
magnetic filed is applied perpendicular to the two dimensional 
electron gas, so that it allows two edge states labeled by 1 
and 2. Electrons are injected from the contact biased at a 
voltage V, red rectangle, and are transmitted to the right- 
hand-side or reflected to the left-hand-side. Two additional 
gate electrodes, represented by blocks (7l and Ur, are placed 
on both sides of the bar and provide a controlled mixing of 
the two channels. Electric current is then measures at the 
exit ports L and R. 



B. Hanbury-Brown and Twiss interferometer 

We shall now demonstrate that it is possible to use 
the witness method to detect the orbital entanglement 
generated in the electronic equivalent of the optical 
Hanbury-Brown and Twiss (HBT) experiment, proposed 
in Ref. [l5|. The electronic HBT interferometer, sketched 
in Fig. [5l consists of a rectangular Hall bar, whose in- 
ner and outer edges are indicated with black lines, and 
four QPCs, labeled with A, B, C and D. Six contacts are 
present, indexed with a number from 1 to 6. Electrons 
are injected from sources 1 and 2, where a voltage V is ap- 
plied against all other contacts which are grounded. Each 
QPC acts as a beam splitter and can be parametrized by 
the following scattering matrix: 



G _( cos^e^ii sin^e*' 
1 ~ I sin 6ie iv ^ cos 9^*™ 



(33) 



At each QPC electrons are transmitted with probabil- 
ity % = cos 2 9i (with i =A,B,C,D) and reflected with 
probability IZi, with IZi = 1 — %. The phases tpj k are 

real and satisfy the condition Lp\ x + ip\ 2 = 77 + P12 + ^21 
in order to guarantee the unitarity of the matrix G;, i.e. 
G\Gi = 1. Along the bar, the electrons follow the edge 
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r 



J 



r 
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FIG. 5: Sketch of the electric Hanbury Brown-Twiss interfer- 
ometer. The edges of the Hall bar are drawn with black lines, 
edges states are represented by red lines and a magnetic flux 
$ threads through the center of the structure. The four QPCs 
are labeled with the letters A, B, C, and D, while the six con- 
tacts are depicted by colored rectangles. Contacts 1 and 2 are 
biased at a voltage V , while all the other are grounded. Elec- 
trons, propagating along the edge states, accumulate phases 
0i, 02, 03 and 4>4- 



states (red lines) in the directions indicated by the ar- 
rows and accumulate phases 0i to 04 (see Fig. They 
can be modulated by varying the length of the paths 
through additional gates. Since there are no interfer- 
ing orbits, no single-particle interference will take place 
and conductance measured at contacts 3 to 6 will be 
phase-insensitive. Because of the indistinguishability of 
the paths taken by the particles injected from sources 1 
and 2, two-particle interference effects arise influencing 
the current cross-correlations measured at the grounded 
contacts. 

Let us briefly recall here the entangled state generated 
in this system by closely following Ref. Il5j. To this aim, it 
is useful to denote the transmission and reflection prob- 
abilities at the QPC C with T c = 1 - Tic = T and at 
the QPC D with T D = 1 - K D = K = 1 - T. The state 
generated by the two independent sources 1 and 2 can 
be written as: |*) = l\ Q<E<eV a[(E)al(E)\0) where |0) 
is the ground state, describing a filled Fermi sea in all 
reservoirs at energies E < and we have omitted the 
spin index. The operator Oj(E) creates an injected elec- 
tron from reservoir j at energy E. In the case of strong 
asymmetry, i.e. 1Z <C 1, it is possible to express the state 
to leading order in \flt as 



|f> = |0> + Vft|*>, 



where 



I*) 



dE 



bl B (E)b 2A (E)-b\ B (E)b 1A (E)] |0> 



(35) 

and b\j (b* 2i ) is the creation operator for electrons coming 
from contact 1 (2) and exiting from the QPCs C and D. 
The second index in the operator indicates towards which 
QPC the electrons are propagating. The state 



|0> 



n 



bl A (E)bl A (E)\0) 



(36) 



0<E<eV 



can be thought as a new vacuum state representing a 
completely filled stream of electrons. The state 
describes an orbitally entangled pair of electron-hole 
excitations^. 

In order to simplify the witness construction it is useful 
to re-write the electron-hole entangled state of Eq. (|3"5|) 
in terms of the vacuum state 10): 



I*) = 



II [&1a(£)4b(£) + 4 A (E)b\ B (E)] |0> 



0<E<eV 

(37) 

Let us now assume that the HBT is subject to dephas- 
ing. Specifically, we consider that the electrons propa- 
gating along the wires besides the phases <pk accumulate 
a random phase Xk, with k = 1, 4. Moreover, this can 
include a randomly fluctuating external magnetic flux $. 
Formally, we introduce the dephasing as follows: 



b\ A 



(E) 

4b(#) 



Te l ^ 1+Xl) b\ K {E) 

Tle l ^ +Xi) b\ B {E) 
7^(02+^) b\ B {E) 



It is quite straightforward then to calculate the density 
matrix p of the outgoing entangled state |\&): 



1 



u(0)eV 



P{E) 



(38) 



0<E<eV 



with 



P(E) 





1 e -i<i>o e -iS o 

e l V^° 1 





(39) 



E 



where we have 



introduced the total phase 
0o = 0i + 02 — 03 — 04 + 27r<f>/<f>o and the dephas- 
ing parameters 5 = x\ + X2 — x% — x\. By assuming 
Gaussian distributions for the random phases, we take 
into account the presence of dephasing by replacing 



e lS with its average value (e 



-i6\ 



in the density 



(34) matrix p. 
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We can then proceed to calculate the witness operator 
w = <S) a< E<eV W e finding that 



We 



10 

-e^° 

-e"^° 

1 



(40) 



whose mean value correctly gives (W) = Tr(pW) = 

2 

— e~ u . We conclude, on the one hand, that the entan- 
gled electron-hole pair is robust against arbitrary degrees 
of dephasing in the HBT and, on the other, that the mea- 
surement of the witness gives an estimate of the degree 
of dephasing afflicting the HBT. 

The witness can be decomposed as follows: 



W E 



1 



it® t) 



■ cos (ft) ((<r x 



'Oy)) 



(41) 



■ sin 0o ((cr x ® <r y ) + (cr y <g> <x x ))] 



The witness can be evaluated by measuring the zero fre- 
quency current cross-correlations s^g(0) of the current 
fluctuations 6I a and 8I@ at contacts a <G 3, 4 and (3 G 5, 6 
(see App. IA 2[) . In fact, although the state exiting the 
QPCs C and D is the one given in Eq. (|34|) , it is easy 
to realize that the only part that contributes to cross- 
correlation measurements at QPCs A and B is indeed 
the entangled electron- hole pair state |^). The vacuum 
state |0), indeed, carries a current from sources to the de- 
tectors but does not contribute to the cross-correlators. 

At zero temperature and in the limit of long measuring 
times r c <C t m -C T tun , the spin correlators appearing in 
Eq. (|4"2")l can be expressed as 



(a-a®b-a) 



2e 3 VTZT 



(4(0) + 4(0) -4(0) 



4» 
(42) 

The expressions for the different values of s°« (0) are re- 



ported in Appendix IA 21 Here we have used the fact 
that in the HBT interferometer the average total cur- 



rents at the QPCs A and B are = 
(I 5 ) = (J 6 ) = ^KV-£. Finally we find 



(h) 



(W) = - cos0o(cos(0o+<5))-sin0o(sui(0o+<5)) 



-TV and 



(43) 



where the last equality is obtained by averaging the vari- 
able S over the Gaussian distribution. 



and triplet states are generated. It has been shown that 
such witness operator can be decomposed into a sum of 
local spin correlators and, in the limit of measuring time 
much smaller than the electron correlation time, its av- 
erage value can be expressed as a function of current 
cross-correlators. This procedure allows the detection of 
entanglement, through the measurement of the witness 
operator, if the latter does not depend on the unknown 
parameters describing the corruption of the state. Being 
a necessary and sufficient criterion for entanglement de- 
tection, the method of the witness operator allowed us 
to prove that spin singlet and triplet electron states are 
robust against the action of pure dephasing in the sys- 
tem. However, entanglement can be destroyed when spin 
relaxation in the conductors is considered. 



For the sake of illustration, we applied the optimal 
witness construction also to two realistic Hall bar sys- 
tems: the quantum point contact tunneling barrier and 
the electronic equivalent of the optical Hanbury-Brown 
and Twiss interferometer, where electron-hole entangled 
pairs are generated. Dephasing was introduced through a 
phenomenological model of randomly distributed phases 
accumulated by the particles propagating through the 
edges channels. Apart from demonstrating the useful- 
ness of the witness method in such setups, we showed 
that the witness detection can be used as a diagnostic 
tool in order to infer information about the "noise" af- 
fecting the system. Unfortunately, a measurable optimal 
witness operator could not be found, using the separabil- 
ity criterion, for the tunneling barrier system in the case 
of arbitrary channel mixing. 



Compared to a Bell inequality test, the detection of 
entanglement by means of optimal witness operators 
presents two remarkable advantages. Firstly, being a suf- 
ficient and necessary criterion, entanglement is always 
detected in the considered class of states. Secondly, in- 
dependently of the degree of dephasing, the evaluation 
of the mean value of the witness requires a fixed set of 
cross-correlated measurements, corresponding to a small 
number (from 10 to 20) of current cross-correlations. A 
drawback is that the method of the witness requires some 
knowledge of the state produced by a given system (in 
terms of an expected class of states) and that the proce- 
dure for constructing a measurable optimal witness op- 
erator is not given a priori. 



V. CONCLUSIONS 

In this paper we discussed how the concept of entangle- 
ment witness can be used in order to detect entanglement 
for electrons in multiterminal mesoscopic conductors af- 
flicted by "noise" . By resorting to the separability cri- 
terion, we constructed an optimal witness operator for a 
prototype setup where maximally entangled spin singlet 



For the future it would be interesting to find other 
schemes for constructing optimal witness operators which 
allow a detection of entanglement even in the cases where 
the separability criterion fails in producing a measurable 
witness. To conclude, we stress that the method of the 
witness operator presented here can be applied to any 
systems and we hope that it can be of practical use in 
experiments. 
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APPENDIX A: CURRENT OPERATOR AND 
CURRENT CROSS-CORRELATORS 

1. Entangler prototype setup 



At zero temperature, the correlators (N^ a N^) can be 
calculated as expectation values over the outgoing state. 
Along the z axis, one finds that 



The current operator for electrons flowing through ter- 
minal A with spin a along the spin quantization axis 
identified by the vector a is defined as^: 



lL(t) = \\ dEdE'e^ E - E '^ H x 
Aa a ( E Ws a {E') - b{ Sa (E)b ASa (E') 



(Al) 



where a}^~ a {E) (b ASa {E)) is the creation operator for in- 
coming (outgoing) electrons with energy E. At zero tem- 
perature, the mean current (J^ a ) can be calculated as the 
expectation value of (|A1[) over the outgoing state, whose 
density matrix is given in Eq. ©. As a result, 



fa) 



e 2 V 



(Nf a ) 



\tA\ 2 \t B f 



(A5) 
(A6) 



where z is the versor of the z axis. Along the x and y 
axes one finds that: 



«U)<(t)> = m^tl + C-ir^cosCxA-XB)] 



(N1 ni) N£ ni) ) 



4 

\tA\ 2 \t B \ 2 



(A2) 

with a = x,y. 



[1 + (-1) S C0S (XA - Xb)] 



where Nj a = b ASa (E)bASa(E) is the number operator. 

The zero-frequency current cross-correlator, defined in 
Eq. (|13[) . can be written in terms of the number operators 
asM: 



4(0) = - {NLwip) 

For the prototype setup of Sec. IIIII we find that 



\t 3 \ 2 



(A3) 



(A4) 



HBT interferometer 



The zero frequency current cross-correlations between 
contacts 3 and 4, and 5 and 6 of the HBT can be cal- 
culated using the parametrization of the QPC given in 
Eq. 



s 36 (0) = 


2e 3 V 


h 


s 45 (0) = 


2e 3 V 


h 


s 35 (0) = 


2e 3 V 


h 


s 46 (0) = 


2e 3 V 


h 



KT 
KT 
KT 
KT 



„A B < 



-A B' 



\9\i9i2) + \9'i<i9n) + 2 I e 

1 A R \ 2 . / A R\2 



i(<j> +5) 


— eT 




A B A B 
3ll52l92lffll 


i(0o+<5) 


— e~ 


i(tpo+S)^ 


A B A B 


i(0o+<5) 


— e~ 


i(4>o+S)^ 


9ll92292l9?2 


i(«fcj+<5) 


— eT 




A B A B 

9i292i9i29n 



(A7) 



where denotes the element of the matrix 

Ga(b) °f Eq. (|33p . Notice, in particular, that such cur- 
rent correlations depend on 9i and the phases <p* fe . 



by setting 



We define the zero frequency cross-correlators srg(0), 

where a,b = x,y,z, a = 3,4 and /3 = 5, 6, by choosing 
different configurations for the QPCs A and B. Indeed, 



1. 9i 

2. 0i 

3. 0i 



0; <pu = 

Ti ^22 



0; f 2 2 = ¥>12 = ¥>21 = 

= 0; </ 22 = tt; <p\ 2 =tp l 21 =0 (A8) 

= J; til = -J; <pii = <Pi2 = 
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we obtain 



1. (.', /. 

2. d = X 

3. G t = Y 



1 

-1 

J_/l 1 
V2 U -1 
J_/l 1 



(A9) 



corresponding, respectively, to the settings z, x and y. 
For example, we find 



4f(o) 



2e 3 V 



TZT 



{X u Y 21 f + {X 2l Y n f 



e 3 F 



-ftT [1 + sin (<f> + 5)] 



(A10) 



X 
Eq 



7 Pjj / 



being the elements of the matrix X (Y) of 
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